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Abstract 

The supersymmetric custodial triplet model adds to the particle content of the MSSM three 
SU{2)l triplet chiral superfields with hypercharge Y = (0,±1). At the superpotential level the 
model respects a global SU{2)l ® SU{2)r symmetry only broken by the Yukawa interactions. The 
pattern of vacuum expectation values of the neutral doublet and triplet scalar fields depends on the 
symmetry pattern of the Higgs soft breaking masses. We study the cases in which this symmetry is 
maintained in the Higgs sector, and in which it is broken only by the two doublets attaining different 
vacuum expectation values. In the former case, the symmetry is spontaneously broken down to the 
vectorial subgroup SU (2)y and the p parameter is protected by the custodial symmetry. However in 
both situations the p parameter is protected at tree level, allowing for light triplet scalars with large 
vacuum expectation values. We find that over a large range of parameter space, a light neutralino 
can supply the correct relic abundance of dark matter either through resonant s-channel triplet 
scalar funnels or well tempering of the Bino with the triplet fermions. Direct detection experiments 
have trouble probing these model points because the custodial symmetry suppresses the coupling 
of the neutralino and the Z and a small Higgsino component of the neutralino suppresses the 
coupling with the Higgs. Likewise the annihilation cross sections for indirect detection lie below 
the Fermi-LAT upper bounds for the different channels. 
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I. INTRODUCTION 


Although the Standard Model (SM) of strong and electroweak interactions is incredibly 
successful in explaining all high- and low-energy particle physics data, it is known to be 
incomplete. One such reason is the astrophysical observation of dark matter leading to the 
belief that there should exist a particle explaining it. However there is no candidate dark 
matter particle in the Standard Model. Moreover cosmic microwave background measure¬ 
ments can be fitted extremely well with a cosmological ACDM model if the relic density of 
dark matter is given by = 0.1187 DEI. 

On the other hand, the recent discovery of the Higgs by the ATLAS [3] and CMS [1] 
experiments seems to point toward a single SU{2)l doublet as being responsible for the 
breaking of the electroweak symmetry. However, this cannot be known for sure without more 
precise measurements of its properties. In fact, beyond the Standard Model proposals predict 
deviations in the couplings of the Higgs compared to the SM values and can alleviate any 
possible future discrepancies between the predicted and observed properties. The simplest 
such models are those with extended Higgs sectors. A drawback for extended Higgs sectors is 
that they can run into trouble with the p parameter if the extra SU (2 )l representations are 
large enoughModels preserving custodial symmetry, for instance, a septet under SU (2 )l 
with hyper charge Y = 2 [MQ] or the Georgi-Machacek (GM) model, which contains a real 
and a complex scalar triplet with Y = (0,1), respectively [5l [8| [TT1 - I36] . can solve this problem 
and keep p = 1 at tree level. However, there is nothing in the previous models protecting 
the scalar masses from large radiative corrections, so all of them suffer from a more severe 
hierarchy problem than the Standard Model. 

Supersymmetry provides a nice solution to both the dark matter problem and the hier¬ 
archy problem. If R parity is assumed, the lightest supersymmetric particle (LSP) is stable. 
For the model to be cosmologically viable, the LSP must be neutral, thus providing a dark 
matter candidate. However, not all dark matter candidates yield the observed relic abun¬ 
dance. Only specific regions of parameter space will allow the lightest neutralino to freeze 
out to the observed relic abundance. In addition to this, the minimal supersymmetric Stan¬ 
dard Model (MSSM) runs into difficulties when trying to ht the observed Higgs mass. One 
method of raising the Higgs mass in supersymmetric models is to extend the Higgs sector 
beyond that of the MSSM. The supersymmetric custodial triplet model (SCTM) [371 EH] 
(a supersymmetric generalization of the GM model) does this by introducing new F-term 
contributions to the tree-level MSSM Higgs mass. Besides, a custodial potential is interest¬ 
ing from a dark matter perspective. The coupling of the Z to the neutralinos vanishes at 
tree level in the custodial limit of the MSSM (for tan/3 = 1), leading to blind spots in the 
spin-dependent dark matter searches [39]. A custodially symmetric extended Higgs sector 
will maintain this property. 

It was argued in Ref. [3B] that a totally custodial situation at the electroweak (EW) scale 
is not favored from a theoretical point of view. The authors examined this issue by imposing 
a global SU{2)l ® SU{2)ji symmetry at a high scale, which would then be broken through 

^ Extra doublets (and singlets) do not suffer from this problem. However representations beyond that do. 
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the renormalization group running of parameters down to the EW scale. Because of the 
influence of the top quark Yukawa coupling the running differentiates the two soft doublet 
masses from each other much more than the three triplet ones among themselves, resulting 
in a much bigger vacuum misalignment in the doublet sector. The difference in the doublet 
sector results in a departure from tan/3 = 1. As the p parameter is affected only by the 
difference in the triplet VEVs, the model still allows the triplets to contribute up to 15% to 
the breaking of the EW symmetry. 

In our study, we take the middle ground between the calculable fully custodial model 
of Ref. [2Z] and the model of Ref. [3H] in which an ultraviolet completion is proposed. We 
assume a Higgs sector with a potential allowing for a noncustodial vacuum, provided that this 
only comes from the ratio of the doublet vacuum expectation values (VEVs), parameterized 
by tan/3. This turns out to be a very good approximation to the situation explored in 
Ref. [38]. Within this approach, we explore the dark matter properties of the model and 
hnd that there are large regions of parameter space where dark matter annihilation in the 
early Universe occurs through the new triplet states. We also examine the direct detection 
consequences of breaking the custodial symmetry along the tan /3 direction and the indirect 
detection bounds. 

The rest of the paper is organized as follows. In Sec. [TT| we briefly review the conditions 
needed to generate the correct relic abundance of neutralino dark matter. Section intro¬ 
duces the model and the benchmark parameters used for our study. The scalar spectrum is 
studied in Sec. W The mixing of the neutralinos is discussed in Sec. |V] which leads into the 
study of the relic abundance of dark matter and direct detection constraints of the model in 
Sec. |Vll We discuss our conclusions in Sec. IVIII 


II. NEUTRALINO DARK MATTER 

In this paper, we focus on dark matter coming from the neutralino sector of the SCTM m 
135] in which the content of the MSSM is extended by three triplets with hypercharge Y = 
(0, ±1). This adds three neutralinos to the four MSSM ones, the Bino, Wino, and the two 
Higgsinos. There are also two new charginos on top of the charged Wino and Higgsino, and 
hnally a doubly charged triplet fermion. We will collectively refer to the fermion components 
of all of the triplet helds as tripletinos. The combination of the neutralinos, charginos, and 
doubly charged tripletinos will be referred to as electroweakinos. 

If the mass parameters of the electroweakinos are well separated, mixing can be neglected 
and the LSP can be a pure gauge eigenstate. The pure Bino does not annihilate enough in 
the early Universe, while both the Wino and Higgsino annihilate easily and need a mass near 
or above a TeV in order to freeze out with the correct relic abundance. If their masses are 
lighter than this, the pure Wino or Higgsino leaves too little dark matter. The pure Wino 
may already be excluded by astrophysical gamma ray searches]^ After constraints from LEP, 
the LHC, and astrophysics are applied, the only pure state that can generate the observed 

^ As SU{2)l triplets, the tripletinos should behave similarly to the Wino in this regard. 
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relic abundance is the Higgsino P0lfl2] . 

To have neutralino dark matter lighter than a TeV and freeze out to the observed relic 
abundance, the LSP must have a large Bino component, and there must be a process which 
helps the LSP to annihilate efficiently in the early universe. There are a few options to 
increase the rate at which the Bino annihilates. 

1. Mixing: If the composition of the LSP contains a substantial amount of Wino, Higgsino 
or Tripletino, the mixing can allow for efficient annihilations. 

2. Coannihilation: Having another supersymmetric particle slightly above the mass of 
the LSP opens the possibility of f-channel annihilations, which can greatly increase 
the annihilation cross section. For there to be enough of the heavier particle around as 
the universe expands and cools down, the mass must not be more than ~ 10% larger 
than the mass of the dark matter (DM) candidate. 

3. Funnel/Resonance: If the mass of the LSP is approximately half the mass of another 
state, the s-channel propagator becomes very large. There is a peak in the annihilation 
cross section, and a corresponding dip in the relic abundance after freeze-out. 

If the LSP is coannihilating with squarks or sleptons, there are strong limits on the model 
from LHC searches. This is due to the production rate of squarks and the relatively clean 
signals for sleptons. In this case, one would expect to hnd the squark or slepton before the 
dark matter candidate. 

In the literature, both mixing and coannihilation among electroweakinos are referred to as 
well tempering |13] . Well tempering implies that there are multiple states around which can 
be produced, which is good for the production cross section of beyond-the-Standard-Model 
states. However, achieving the correct relic abundance requires the splitting to be small, 
which makes detection difficult. There have recently been studies on detecting electroweaki¬ 
nos with small splittings at colliders |11H52]. 

The resonant/funnel annihilations of the LSP do not need extra particles at the same 
mass, but instead at nearly twice the mass of the dark matter particle. In the MSSM, the 
funnel particle can be either of the (PP-even Higgs, or CP-odd Higgs (H°) |53H59]. 

The dark matter particle itself cannot be detected at colliders, which implies the way to 
look for such a model is through the heavier states. However, searches for neutral scalars 
are difficult, as exemplihed by the long search for the Higgs. As will be shown later, in 
the SCTM, the triplet scalars provide a resonant channel over much of the parameter space. 
Because of to the degeneracy of states in the custodial situation EH. there are charged states 
near the neutral funnel that could aid in discovery. 


III. MODEL 

As in Ref. 133 , we will construct the supersymmetric Higgs sector manifestly invariant 
under SU{2)l ® SU{2)r. The MSSM Higgs sector Hi and H 2 with respective hypercharges 
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F = (-1/2, 1/2) , 


Hi = 


HI 


Ho = 


HI 


is complemented with SU{2)l triplets, Sy, with hypercharges Y = (—1, 0, 1) 


S-i = 



X 


X 

x/2 


Sn= 



Si = 


_ V2 






( 1 ) 

( 2 ) 


where Q = T^l + Y- 

The two doublets and the three triplets are organized under SU{2)l ® SU{2)r as H = 
( 2 , 2 ), and A = (3, 3) where 


H = 




(3) 


and T^ji = Y. The invariant products for doublets A ■ B = A°'eabB^ and antidoublets 
A - B = Aae^’^Bc are defined by 621 = = 1 . 

The SU{2)l ® SU{2)r invariant superpotential is defined as 

Wo = \H ■ hH + + -H ■ H + (4) 

3 2 2 


and the total potential 


where 


V = Vf + Vd + F,oft 


(5) 


Koft = tr |Sip + tr |S_ip + tr |SoP 

+ I ^m\H ■ H + ^i?A tr + AxH ■ AH + tr A^ + h.c. | (6) 

Note that the potential we just wrote is the same as in Ref .1371 but with noncustodial 
soft masses. They will be used to satisfy the equations of motion. The neutral components 
of all the helds can be parametrized by 

X = ^{vx + Xr + iXi), X = ( 7 ) 

By imposing 

ui = x/^cos/duH, V 2 = \/2sm. (5 vh and v^p = = v^p = va, (8) 

the custodial symmetry is only broken in the vacuum by tan 13. For the rest of the paper we 
refer to tan (3 = 1 as the custodial case and tan (3 ^ 1 as the noncustodial case. To set the 
Z mass, the total VEV must be 

= (246 GeV)^ = 2u| + 8^. (9) 
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The p parameter is not affected if custodial symmetry is broken in this way. This will generate 
hve equations of minimum that we will use to solve for the values of the soft masses. The 
minimization conditions can be found in Appendix While Eq. ( |A 6 [ ) (and the previous 
noncustodial equations) provide the necessary conditions for an extremum of the vacuum, 
they do not guarantee a minimum. To guarantee the appearance of a minimum, it is a 
sufficient condition that the determinant of the Hessian at the origin be negative. In the 
custodial case, to leading order in small ua, this condition can be expressed as 

3 

A(2/i — Pa) — Aa > 0 and ~ < 0. (10) 

The left equation sets the allowed relative size between the doublet and triplet supersymmet¬ 
ric masses. The right equation will be very important as we scan across the parameter space. 
For a fixed value of rUg, the right equation sets a maximum value for The A parameter 

will be used to raise the tree-level Higgs mass, so there exist regions where the Higgs mass 
cannot be achieved with light stops while keeping the potential correctly minimized (i.e. not 
getting tachyonic states). 

To begin a study of the dark matter properties of the model, we first choose a set of 
benchmark values, given by 


Ag = 0.35, 
m 3 = 500 GeV, 

Ba = -(500 GeV)^ (11) 

Ax = Ax^ = At = Ai) = A-r = 0 , 

uIq^ = 800 GeV, and = 700 GeV, 

where other scalar soft masses have been decoupled and the ones corresponding to Higgs 
multiplets are determined by the minimization conditions. The SGTM triplet F terms yield 
a large tree-level Higgs mass, so smaller one-loop corrections are needed. This is the reason 
for our choice of relatively light stops albeit above the current experimental limits. The value 
of A 3 will not have much of an effect. We are considering the case m| = |i?A| for simplicity. 
Values of m 3 and Ha around those in Eq. (0 should provide similar results. Larger values 
for m 3 or Ha will decouple the heavy scalars more, and in addition will affect how large ua 
can be in the minimization of the potential. Similarly, we choose to examine the case in 
which all of the trilinear terms are zero to help ensure that the EW vacuum is the deepest 
one. This leaves p, ps. A, and ua as the remaining free parameters to study. 


IV. SCALAR MASSES 

There is a total of hve GH-even, hve GH-odd, six singly charged, and two doubly charged 
Higgs scalar helds in this model. The mass matrices for all of these states are cumbersome, 
and not entirely enlightening. In Ref. |37], the helds were decomposed into the SU{2)v 
custodial basis, which is also directly related to the mass basis for tan/? = 1 , up to small 
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hypercharge breaking effects. This notation can be helpful for showing how many charged 
states will be in the proximity of a neutral state. However, since we will be examining 
both the custodial and the noncustodial setups of the model, we will not use this notation. 
Instead, we will work with mass eigenstates. After removing the Goldstone bosons, they will 
be denoted as A? 4 , and T^^ 2 - 

To study the dark matter annihilation in the model, we are really only interested in the 
spectrum of the lightest neutral scalars rather than the charged components. Annihilating 
the neutralino through a resonance of the Higgs or the heavy Higgs has been shown before 
in the MSSM. As a new feature of this model, there are substantial regions of parameter 
space in the SCTM where the annihilation can proceed through a tripletlike resonance. To 
do this, the soft masses of the triplets must not be too large. Upon close examination of the 
minimization conditions for and in Eqs. (A5), (A3), and (A4) respectively, 

we see that there is a piece that scales as vjj/vA for each soft mass. Smaller values for va 
yield large soft masses for the triplets, decreasing the chance of annihilating through the 
triplet funnel. 

In the decoupling limit of the MSSM, the tree-level Higgs mass goes as cos^(2/3). 

Because of this, for tan/? = 1, there is no tree-level contribution to the mass of the Higgs 
from the MSSM parts of our model. Instead, the mass at tree level in the decoupling limit 
comes only from the triplet F terms, and is given (at leading order in va) by 


m. 


= 3X^v 


‘'h\tanf3=i- ( 12 ) 

We also examine the model in which tan /? 7 ^ 1 but is still small. The tree-level Higgs mass can 
no longer be written in a simple form. However, we comment that there are now MSSM con¬ 
tributions to the mass, and the triplet F terms contribute as (4 cos^ /? 4 - 4sin^ /? -|- sin^ 2/3). 
The SCTM allows for large tree-level contributions to the Higgs mass with no need of large 
one-loop corrections, and thus no need for heavy stops. 

The dominant radiative corrections to the Higgs mass depend on the top Yukawa coupling, 
defined as 

, rrit rrit rrit , , 

ht = —^ J ^ J (13) 


V2/\Pi 


sin/3 uh sin/3y'(r;2 - 8 n|) /2 
They have been proven to be sizeable in the context of the MSSM, without jeopardizing 
perturbativity, as the top Yukawa coupling does not enter the Higgs mass at the tree level. 
In fact in the SCTM, increasing va increases the top Yukawa, which increases radiative 
corrections to the Higgs mass. In our study, we take the dominant one-loop corrections 
found in Ref. [60], we use 700 CeV for the right-handed soft mass and 800 CeV for the 
left-handed soft mass. These were chosen to be slightly above the current experimental 
bounds, regardless of the mass of the lightest neutralino. Raising the masses of the stops 
will not affect our dark matter results, only worsen the fine-tuning of the model. Note that 
even though the stop masses and At are fixed in the study, changing fi and tan /? affects the 
mixing and thus the one-loop contributions to the Higgs mass. 


^ We will neglect radiative corrections proportional to as the parameter A affects the Higgs mass at the 
tree level and thus the corresponding radiative corrections are constrained to be small by perturbativity. 
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Figure 1. Top row: Maximal values of ua that allow A to set the Higgs mass to 125 GeV and yield 
a minimized potential as a function of Bottom: Value of A needed to attain the observed Higgs 
mass for ua = 10 GeV. The triplet supersymmetric mass is set to = 250 GeV, and the other 
values are as in Eq. ©• 


By fixing the mass of the stops, the only way to alter the mass of the Higgs is through 
the remaining parameters, /i, /is, A, and ua- To study the effect of the triplet states on dark 
matter, we examine the case in which either the doublet- or the tripletlike fermions are 
lighter. We £x /ia = 250 GeV and scan over the values of /x. 

Recall that in order to achieve a minimum of the potential from the minimization condi¬ 


tions in Eqs. (A1 )-(A5), rather than a saddle point, there exist constraints on the relationship 


between /i and /xa beyond Eq. (10). However, the latter equation might give us some in¬ 


tuition on which /x and /xa values we can take since a saddle point at the origin forces the 
potential to have a minimum. When tan /3 = 1 and A;, = 0 the equation simplihes to 
2 /x > /Xa and we see that we cannot look at regions where /xa is signihcantly heavier than /x 
and still minimize the potential]^ 


Of course, this does not mean that triplets cannot be decoupled supersymmetrically. The limit /xa —t oo 
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Once yU and /ia are fixed, we have to chose A and va- We do this in two different ways: 

1. To maximize the value of va, we start with va = 30 GeV, which we take as the upper 
limit as suggested by the analysis of Ref. [SB]- We then scan over A to set the Higgs 
mass (including radiative corrections). Once the lightest GP-even Higgs has a mass of 
125 GeV, we examine the rest of the spectrum. If other scalars have gone tachyonic, 
or the value of A needed is greater than 0.75, this value of va is excluded. We then 
lower Va and repeat the process until a 125 GeV Higgs is obtained and the vacuum 
minimized. The resulting values of va and A are plotted in the top row of Fig. [^over 
a range of /i. 

2. The other option is to keep the value of va constant as we scan across fi. The region 
of /i that can yield the correct Higgs mass and successfully minimize the potential is 
smaller for large values of va- Because of this, we set va = 10 GeV for our study 
of this method. The lower panel of Fig. [T] displays the values of A needed for both 
tan/3 = 1 and tan/? = 2. Note that tan/? = 2 needs smaller values of A because there 
are tree-level MSSM contributions to the Higgs mass. This allows for a larger range of 
fi than the tan /? = 1 case. 

The spectrum of the light neutral scalars is plotted in Fig. for va = 10 GeV and va 
maximized in the left and right panels respectively. When tan/? = 1, shown in the upper 
panels, the scalars H 2 , H?, and have similar masses, which increase as a function of fi. 
The lightest that these scalars can be is ~ 300 GeV. The other neutral scalars all have 
masses greater than 600 GeV and therefore are not shown in the plots. In the lower panels, 
the same spectra are shown for tan /? = 2. In this case, both 77° are nearly degenerate 

in mass, and much lower in mass than when tan /? = 1. This partially comes from the smaller 
value of A needed to raise the Higgs mass for tan/? = 2. Gonversely, the mass of 77° does 
not change much between the two choices of tan /?. If the maximum va is chosen instead of 
using the constant va = 10 GeV, the masses of 77°, H?, and 77° will drop. The separation 
of the states will also depend on va so increasing it helps to remove the degeneracy of the 
scalars. 

We do not perform any collider constraints on searches for these extra possible scalars. 
However, we see that the model allows for some to be very light. A dedicated search could 
therefore exclude large regions of parameter space in a quicker and more conclusive way than 
either Higgs precision measurements or direct detection experiments. 


V. NEUTRALINO MIXING 


The addition of three triplet chiral superhelds adds to the neutralino content of the model 
three extra states. The mass Lagrangian in the basis ^jJ^ = ^.B°, 1V°, 77°, 77°, 0°, ;y°,-0°^ is 
then 


-^neutralino mass — 2^^ ^ IVI'0 -|-C.C., 


( 14 ) 


yields the MSSM, in which case Eq. (101 does not apply. 
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Figure 2. The left panels show the spectrum of the neutral light scalars when ua = 10 GeV and A 
is changed to set the Higgs mass. The right panels use the maximum allowed value for ua for each 
/i value. The upper (lower) panels contain tan/3 = 1 (tan/3 = 2). Changing tan/3 greatly affects 
the masses of and A^, but mass is similar for both choices. 


where 


M = 


/ 

Ml 

0 


^g'sfjVH 

0 

-g'vA 

g'vA 


0 

M 2 

^g2CpVH 

-j^g2S0VH 

0 

g2VA 

-g2VA 


-j^g'yvH 

7^92CpVH 

— y/2XvA 

-j^XvA-fi 

-XspVH 

0 

—2XCi3Vh 
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Figure 3. Composition of the LSP in terms of gauge eigenstates. The top row shows tan/3 = 1 and 
the bottom shows tan/3 = 2. The columns correspond to /r = (200,250,400) GeV, respectively, 
while the tripletino mass is set to /ta = 250 GeV. The Wino has been decoupled with M 2 = 1 TeV. 
Note in the top middle and top right plots the presence of a tripletlike eigenvalue, which is totally 
decoupled from the rest of the neutralino mass matrix, made out of only ijj, (p, and y. It corresponds 
to an SU{2)v 5-plet in the custodial basis. 


and Sjs and are shorthand for the sine and cosine of 13, respectively. 

Overall, the masses are controlled by Mi, M 2 , 11 , and /ia for the Bino, Wino, Higgsinos, 
and tripletinos, respectively. There are additional contribntions to the masses and mixings 
scaling with either vh or va- To provide a good dark matter candidate, we want the LSP 
to be the lightest nentralino; its composition will then determine the annihilation and direct 
detection cross sections. 

The composition of the LSP in terms of the gange eigenstates is shown in Fig. for 
the case in which the VEV of the triplets is constant 10 GeV and tan/3 = 1(2) in the top 
(bottom) row. The left panels have the Higgsino-like states lighter than the tripletino ones, 
nsing /i = 200 GeV and /ia = 250 GeV. The middle panel has both the Higgsino and 
tripletino masses set to p = /ia = 250 GeV. Finally, the right panel examines when the 
triplet states are lighter than the Higgsino, with fi = 400 GeV and /ia still at 250 GeV. To 
simplify the sitnation as much as possible, we decouple the Wino by setting M 2 = 1 TeV. 

In the custodial situation, the doublet components of the LSP are equal and the triplet 
components are separately equal over most of the parameter space. The tan [3 = 2 case has 
each Higgsino and tripletino contributing differently to the LSP. Despite the complexity of 
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the plots, there are a few overarching trends. 

In Sec. [TT[ we argued that the Bino component of the LSP must dominate in order to 
achieve the correct relic abundance of dark matter. The interesting regions to examine 
in the compositions plots are then Mi < p,pA- In this region, even when p > /ia, the 
second-largest component of the LSP is Higgsino rather than tripletino which is true even 
for quite large values of the Higgsino mass. This is due to the mixing of the Bino with 
the Higgsinos or tripletinos, which comes from off-diagonal terms weighted with vh or va, 
respectively. Because of Eq. S> ua, and the Higgsino mass needs to be much larger 

than the tripletino mass in order for the triplet contribution to the LSP to be larger than 
the Higgsino component. So even though the mass of the Higgsino can be larger than the 
tripletino mass, the mixing of the Bino with the Higgsino can be what causes the correct 
annihilation rate. 

As /i is further increased, the amount of Higgsino in the LSP drops past the point where 
mixing alone can yield the correct relic abundance. Looking only at regions where Mi < /ia, 
we see that the triplet states do not contribute much to the LSP. By removing the Higgsino, 
the LSP is made more pure Bino, rather than increasing the triplet amount. The only 
possibility of well tempering for this will then require coannihilations of the Bino-like LSP 
with a tripletlike state. 


VI. DARK MATTER 

To examine the dark matter of the SCTM the model was implemented into SARAH 
[SIHS5]- With this, a code was generated for SPheno [BSl EZ] and CalcHep pS]- The 
SPheno code calculates the spectrum, outputting a parameter card that can be read by 
MicrOMEGAs 3 [On]. The program MicrOMEGAs 3 uses the GalcHep code to cal¬ 
culate the dark matter properties. 


A. Thermal relic density 


For each of the choices of tan (3 and the method of picking va, we scan over the possible 
fi values for /ia = 250 GeV, using 50 GeV step sizes. At each point in /i, we then scan over 
Ml to hnd the Bino masses that yield the correct relic abundance of dark matter. We start 
with Ml = 40 GeV and take 1 GeV steps until Mi > 100 GeV, at which point a 5 GeV step 
is used to save on computing time. 

Figure shows an example of the relic abundance calculated at each Mi value for the 
point /i = 200 GeV, va = 10, and tan/3 = 1. The gray line marks VLh? = 0.1187, the 
observed relic abundance in the Universe PI. The scalar masses do not depend on the Mi 
value and are given by 


mijo = 125 GeV, = 299 GeV, 

rrijj^o = 325 GeV, rnffo = 337 GeV, and others > 700 GeV. 


(16) 
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Relic abundance 
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Figure 4. Relic abundance for the model with /r = 200 GeV, /ta = 250 GeV, ua = 10 GeV, and 
tan/? = 1. The gray line marks the observed relic abundance in the Universe today. As the mass of 
the LSP crosses over half the mass of one of the scalars in the model, the annihilation cross section 
greatly increases, leading to lower relic abundances. When the LSP mass gets close to the mass 
of the Higgsino, the mixing and coannihilations take over and the relic abundance stays below the 
observed value. 


Three dips in the relic abundance are seen in the plot corresponding to the funnel, the 
H 2 funnel, and one for the nearly degenerate and states occurring when the Bino 
mass is roughly half the scalar mass. There are three Mi values of this model point that 
yield the correct relic abundance. The first two correspond to going into and out of the 
lightest Higgs funnel, and the third one is at the start of the iL® funnel. However, the next 
funnels corresponding to A\ and are close together, so the effect of having multiple nearly 
resonant s-channel annihilations keeps the relic abundance below the observed value. This 
runs into the region where Mi > /i and the Higgsino becomes the LSP, leaving not enough 
dark matter in the current Universe. 

For each p value in our model scans, we do the same process. Whenever the relic abun¬ 
dance at one Mi value crosses from one side of the observed value to the other at the next 
Ml step, we do a more dedicated scan to find the Mi value to a higher degree of accuracy. 
We then classify the point according to the process that is driving the annihilations by com¬ 
paring the LSP mass to half the mass of the scalars or 10% higher than the LSP mass with 
that of the next-lightest electroweakino. The piece giving the minimum of 
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(17) 

yields a classihcation of the given scalar funnel or well tempering. This classihcation is only 
an approximation of what is actually causing the annihilations. In the nonrelativistic limit. 
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annihilations through scalars occur through the p-wave, while pseudoscalars occur through 
the s-wave. Thus, when is close in mass to either if® -^ 3 ) the classihcation scheme could 
point to the scalar instead of the pseudoscalar, even though the pseudoscalar contribution 
is larger. In addition, when the funnels are close to the well tempered region both process 
can be responsible for the annihilation. 

The results of the classihcations are plotted in Fig. [^for the different model choices in 
the m^o vs. /i plane. The LSP is mostly Bino, so Mi ~ m^o. The triplet scalars can be very 
light for tan /3 = 2 or if tan /3 = 1 when the VEV of the triplets takes on the maximum value 
allowed. Recall that Fig. shows that these masses increase as a function of /r. As such, 
the funnels for the Tripletlike, H 2 , A^, and H^, scalars smoothly transition up to the point 
where well tempering happens at a lighter mass than needed for a triplet funnel. 

For every model choice examined, there is an Mi value that will yield the correct relic 
abundance either through a Tripletlike scalar or well tempering. When /i is large enough 
that the triplet scalars funnels are not possible, the Higgsinos are heavy enough that the 
well tempering is not caused by Bino-Higgsino mixing but instead by coannihilations with 
the triplet fermions. Thus each model point examined is capable of setting the correct relic 
abundance using particles beyond the MSSM content. 

The large VEV of the triplets allows for the triplet scalars to be light. The lightness of 
these scalars is what allows the model points examined to always be able to set the relic 
abundance using either the triplet scalar funnels or the triplet fermions. However, lowering 
the triplet VEV, ua, raising the triplet supersymmetric mass, /xa, or lowering the Wino 
mass, M 2 , can disturb the possibility of achieving the correct relic abundance through a 
triplet state. The MSSM limit of the model takes the VEV of the triplets to zero. In this 
case, the triplet scalar soft masses go to inhnity and do not contribute to the annihilations]^ 
The Higgsino alone satisfies the correct relic abundance if its mass is ~ 1.1 TeV. As such, if 
/xa is much larger than that, the triplet fermions cannot play a role in well tempering. Such 
a large value of fi would also keep the triplet scalars heavy, so such a case would have no way 
of using the triplet superheld to set the relic abundance. Finally, the Wino has been raised 
above the mass of the Higgsinos and tripletinos for this study. Bino-Wino well tempering 
can also be done if Mi ~ M 2 < /i,/iA- In this case the relic abundance could be set before 
the triplets have a chance to affect things. 


B. Direct detection 

There have been many experimental searches for the direct detection of dark matter. For 
the mass ranges considered here, the Particle Data Group cni shows that the best limits 
are currently coming from the LUX Collaboration for spin-independent searches and the 
COUPP Collaboration [73l [71] for spin-dependent measurements. Super-Kamiokande na 
and IceCube Hg EZI have better spin-dependent exclusions, but are indirect constraints 
that rely on the annihilation of dark matter in the current Universe and depend on the the 

® This also happens in triplet models in which the p parameter is not protected by a custodial symmetry, as 
the triplet extension of the MSSM E], and ua is strongly constrained by electroweak precision observables. 


14 



\/a=10 GeV; tan;S=1 



\/A=max; tan;S=1 



250 


\/a=10 GeV; tan/3=2 


\/A=max; tan/3=2 



200 


400 600 

MGeV] 


800 1000 



^J[GeV] 


• Funnel 

■ Funnel 

♦ Funnel 
▲ Fl§ Funnel 

T Well Tempering 


Figure 5. Points that yield the correct relic abundance of dark matter. The upper row is for the 
custodial case, while the lower has tan (5 = 2. The left panels keep ua constant, and the right panels 
use the maximum allowed value for ua for each /r value. The points are labelled corresponding to 
which annihilation channel dominates in the early Universe. 


byproducts of the annihilation that change as the LSP composition changes. We then only 
compare our results with the LUX and COUPP constraints. 

The spin-independent cross sections for the points satisfying the correct relic abundance 
are shown in Fig. The micrOMEGAs 3 output provides both the cross section of the 
dark matter with a proton and a neutron; we take the maximum of these. The points are 
marked in the same fashion as Fig. to show how the relic abundance is being achieved. The 
upper (lower) panels show tan /? = 1 (2) while the left and right panels display ua = 10 GeV 
and when ua is maximized at each point, respectively. The shaded blue region is excluded 
by the LUX bound, and the dashed blue line is the projected sensitivity of LUX. 

The spin-independent cross section is mediated by the doublet scalars. There is not much 
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Figure 6. Spin-independent dark matter nucleon cross sections. Each point meets the correct relic 
abundance with the annihilation mode marked. Points with smaller Higgsino components have a 
lower spin-independent cross section. 


difference between the tan /S = 1 and tan (3 = 2 models in terms of the cross sections. For 
Vi\ = 10 GeV, both have a region where the dark matter mass is between 100 and 200 GeV 
which can be excluded by LUX. The points are achieved through a triplet funnel, and to get 
masses in this range for the LSP, the values of /i are low. Referring back to Fig. low values 
for /i and Mi give the LSP a moderate Higgsino component. This Higgsino component is 
what drives the nuclear cross sections to be so large. The cross sections are lower when the 
maximum value of va is used. In this case, there are few points that are currently excluded 
by LUX. The larger value of va lowers the masses of the Tripletlike scalars. This pushes 
the triplet funnels and the well-tempering regions to larger values of fi, further decreasing 
the Higgsino component and the spin-independent cross section. Fortunately, there are still 
many points that can be probed by LUX in the future. However, the points which are well 
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Figure 7. Spin-dependent dark matter nucleon cross sections. Each point meets the correct relic 
abundance with the annihilation mode marked. The parity-violating Z couplings vanish in the 
custodial case. 


tempered through Bino-tripletino coannihilations remain under the projected bound, due to 
the minimal Higgsino component of the LSP. 

The spin-dependent interactions are mediated by the Z boson and the cross sections are 
shown in Fig. The panels use the same labelling as Figs. |^and[^ In the custodial case, with 
tan /9 = 1, the mass eigenstates of both the fermions and the scalars of the Higgs doublet and 
triplet superhelds form representations of SU{2)v- The parity-violating Z coupling therefore 
vanishes in this case. And while this is also true in the MSSM for tan/? = 1, the SCTM 
provides motivation for this choice of tan /?. The model points examined for tan /? = 2 
no longer have vanishing Z couplings with the LSP. The cross sections are much larger in 
this case, particularly for the well-tempered points, which have low spin-independent cross 
sections. However, even these large cross sections are still ~ 2 orders of magnitude below 
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the COUPP bound. 


C. Indirect detection 


The direct detection experiments rely on dark matter interacting with detectors on Earth. 
It is also possible to observe astrophysical objects containing large dark matter densities. In 
these regions of space, the LSP can still annihilate. The annihilation does not occur through 
a diphoton process, which would lead to a monochromatic signal. Instead, experiments must 
search for photons coming from the byproducts of the annihilation. 

The annihilation cross section in the current Universe can be much different than in the 
early Universe. Scalar funnels (not pseudo) are velocity suppressed in the nonrelativistic 
limit. As the temperature has cooled since freeze-out, the annihilations proceeding through 
scalars should be signihcantly smaller than the ~ 3 x 10“^®cm^/sec needed at freeze-out. For 
the well tempering through coannihilations, the coannihilating particle is no longer around 
in the current Universe, so we expect the annihilation cross section to be lower now as well. 

The Fermi-LAT Collaboration [781179] has placed limits on the annihilation cross section of 
dark matter from the observation of satellite galaxies. The limits are framed in the context of 
the annihilations proceeding 100% of the time through either the e’''e“, , uu, bb, 

or W~^W~ channel. In Fig. the resulting limits are plotted with our model points yielding 
the correct relic abundance. The upper (lower) panels show tan/? = 1 (2), while the left and 
right panels display ua = 10 GeV and when va is maximized at each point, respectively. A 
few points for the tan (3 = 2 case are possibly excluded by these searches. However, these 
each have the largest annihilation channel being x^Xi —>■ H^Z. The spectrum of photons 
coming from the decays of the and Z will not map directly onto any of the Fermi-LAT 
limits. The fact that the SCTM has more neutral Higgs funnels opens this possibility of 
having different annihilation modes. A more detailed study would therefore be needed in 
order to conclusively exclude points from the SCTM due to indirect constraints. 

We also note that some of the points marked as annihilating through the pseudoscalar 

have particularly large annihilations in the current Universe. These interesting points 
have A? very close in mass to either and there are interference effects in the early 

universe keeping the annihilation cross section small enough. In the current Universe, when 
the scalars do not play as much of a role, the annihilations proceed with less interference. 
Similarly, many points marked as funnels seem to have annihilation rates larger than 
expected in the current Universe. If the rates are scaled up by the larger velocity at freeze- 
out, the annihilation rate would seem to be too large. However, these points he close to the 
well-tempered region, so it is likely that a simple classihcation does not work well for points 
where both processes are important. 
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Figure 8. Annihilation cross section times velocity of dark matter in the Galaxy in the current 
Universe. Each point meets the correct relic abundance with the annihilation mode in the early 
Universe marked. The lines mark the limits assuming the annihilation occurs 100% of the time 
through the given channel, each resulting in different spectra of photons measured here on Earth. 


VII. DISCUSSION AND CONCLUSIONS 

We have studied a supersymmetric model in which the Higgs sector of the superpotential 
is extended by three SU{2)l triplet fields and is manifestly invariant under SU{2)SU{2)r. 
When the triplet VEVs are aligned, the custodial symmetry of this setup allows the triplet 
helds to be light and develop large vacuum expectation values with no effect on the p pa¬ 
rameter but participating in the electroweak symmetry breaking. Therefore, we only allow 
this custodial setup to be broken by tan /?, the ratio of the doublet VEVs. This has an effect 
on the dark matter phenomenology, but does not alter the electroweak precision constraints. 
The F terms coming from the triplet helds help to raise the mass of the Higgs to its observed 
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value, lowering the needed mass of the stops. 

We studied the case in which the lightest supersymmetric particle is a neutralino, which 
with R parity gives a stable dark matter candidate. For the dark matter candidate to yield 
the correct relic abundance, it must have a large Bino component to not annihilate too 
quickly in the early Universe. Well tempering mixes the Bino with either the Higgsino, 
Wino, or tripletino in just the right amount to give the observed relic abundance. If the 
Bino component is too large, dark matter does not annihilate quickly enough in the early 
Universe, unless the mass of the dark matter particle is about half the mass of a boson. We 
fonnd that the triplet scalars or the triplet fermions can play a role in the annihilation of dark 
matter in the early Universe over a large range of values for the Higgsino mass parameter 
/i. We compared the model points giving the correct relic abundance with the cnrrent best 
direct detection limits. The points with low /i valnes have at least a moderate Higgsino 
component and have either been exclnded already or can be discovered in fnture results. At 
large values of /i, the light triplet states still provide an efficient means of annihilating the 
dark matter, bnt hope of a direct detection of dark matter is lost. This is motivation for a 
detailed study of the LHC phenomenology of these models. 

For the stndy of dark matter, we were only concerned with the nentral triplet states. 
However, the triplets contain charged states. In fact, the second-lightest CP-even Higgs, 
H 2 , is close in mass to both the lightest charged and donbly charged scalars. In addition, 
the mixing of the Tripletlike fermions leaves charged and doubly charged states very near in 
mass to the lightest nentral one. If the relic abundance of dark matter relies on these light 
states, they shonld be accessible at the LHC. While a detailed study is beyond the scope 
of this paper, a dedicated study of methods for searching for donbly charged fermions and 
scalars could offer valuable constraints on models such as these with exotic particle content. 

To conclude, the SCTM helps raise the mass of the Higgs through extra F terms. There 
are large regions where the lightest Higgs is consistent with cnrrent observations, and able to 
explain any deviations in fnture measurements. It also offers new methods to annihilate dark 
matter in the early Universe throngh triplet fermion coannihilations or triplet scalar portals. 
The triplet fermions have weak detector bounds, if any, as compared to coannihilations with 
squarks or sleptons. The triplet scalar funnels can also be light, and have charged partners 
making them easier to search for than plain MSSM nentral Higgs funnels. Only a dedicated 
LHC search for the triplet fermions or scalars could offer the most conclusive constraints on 
the model. 
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Appendix A: Minimization conditions 
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